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Abstract 

A proper discretization of the logistic differential equation, which is 
preserving these two distinct equilibrium solutions and their unstability 
and stability, suggest that we need to examine the time delay of the logistic 
map. According to Murray [2], the effect of delay in models is "usually" 
to increase the potential for instability However the word "usually" is 
really ambiguous. In this paper, we mathematically formulate and prove 
the two conjectures about stability and instability. 



1 Introduction 

The logistic differential equation (cf. [3], etc . . . ) is given by 



dec f X . 

rxil- — ), (1.1) 



dt \ K 

where r and K stand for the reproduction rate and the carrying capacity for the 
population, respectively. Thus, let us assume that r is a positive parameter and K 
is any fixed positive constant. Here, we note that the equilibrium solution Xi = 
is unstable for any r > 0, and the equilibrium solution X2 = K is stable for any 
r > 0. 

Let us find a proper discretization of (1.1), which preserves these two distinct 
equilibrium solutions and their instability and stability. First, we consider the 
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following discretization: 



x(t + At)-x(t) x(t)\ 



Xn+l — Xn ( X, 

rXn 1 - 



X (1.2) 

h '""V^ K 

(forward difference) 

where X^+i := x{t + At), X„ := x{t), and h:= At{>0). 
This becomes a first order nonlinear difference equation: 

X„+i = (r/i + l)X„-— X^. (1.3) 

We note that the fixed point X = is unstable for any r > 0, and the fixed point 
X = is stable if < r < |, and unstable if | < r. Therefore, though the 
non-trivial equilibrium solution of (1.1) is stable for any r > 0, the non-trivial 
fixed point of (1.3) is stable only for < r < |. In this sense, (1.3) is not a proper 
discretization of (1.1). 

Next, we consider the following discretization: 

x{t + M)-x{t) x{t + At)\ 
^ ^ V y _ ^^(^^-^ ^ — )_ as At ^ 



^ Xn+l — Xn ^ f X„+i ^ ( • ) 

^ h = ^^"i^-^ 

(forward difference) 

where := x{t + At), X„ := x{t), and h := At(> 0). 
This becomes a first order nonlinear difference equation: 

X„,. = ^1^. (1.5) 

We note that the fixed point X = is unstable for any r > 0, and the fixed point 
X = K h stable any r > 0. Then, the non-trivial fixed point of (1.5) is stable for 
any r > 0. In this sense, we see that (1.5) is a proper discretization of (1.1). 

Observe that if the term X^+i on the right-hand side of (1.4) is replaced by 
Xn, then the range of stability of the trivial fixed point does not change and the 
range of stability of the non-trivial fixed point shrinks from (0, +cxo) to (0, |). 

If we extend this observation, we can naturally guess that if the term Xn^i on 
the right-hand side of (1.4) is replaced by X„_^, then the range of stability of the 
trivial fixed point does not change and the range of stability of the non-trivial fixed 
point shrinks further as r increases. 
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According to Murray [2], the effect of delay in models is "usually" to increase 
the potential for instability. However the word "usually" is really ambiguous. 

In this paper, we mathematically formulate and prove the following conjec- 
tures: for the discrete logistic equation which includes the time delay term X^-t, 
(i) the range of stability of the trivial fixed point does not change as we change the 
value of r; and (ii) the range of stability of the non-trivial fixed point constricts 
gradually as we increase the value of r. 



2 Discrete Logistic Equation with Time Delay 



Equation (1.1) is not realistic in modeling a population with distinct breeding 
seasons for it incorporates instantaneous response to changes in the environment. 
In reality, this response to changes usually takes effect after some time delay or 
time lag. 

In making things more realistic in modeling populations, we shall consider 

another equation which will be the main equation of this study. This equation is 
similar to the discrete logistic equation, however, it incorporates the time delay 
between increased deaths and the resulting decrease in population reproduction. 
This equation is known to be the Discrete Logistic Equation with Time Delay [3], 
which is given by 



Again, r and K are the reproduction rate of the population and carrying capacity 
respectively. We denote the time delay to be Equation (2.1), similar to Equation 
(1.1), has two distinct fixed points, namely the trivial fixed point denoted by Xi = 
(0, 0, ■ ■ ■ , 0), and the non-trivial fixed point given by X2 = {K, K, - ■ ■ , K) . 



Futhermore, we can guess that the stability of the population size of equation 
(2.1) also depends on the value of r as equation (1.1). Equation (2.1) is applied 
by May on A.J. Nicholson's experimental data concerning Australian blowflies 
[3], a pest in the Australian sheep industry. May well-approximated the ex- 
perimental data Nicholson had gathered using equation (2.1) and is shown by 
Figure 3.8.1 in [3]. The parameter values used for the output of the model are 
r = Q.lQ6days-\ K = 2.8 x 10^ flies and r = 17 days. 




(2.1) 



r+1 



r+1 



3 Stability of the trivial fixed point 



In this chapter, we summarizes the results we obtained from the computations. 
From the Table 1 below, notice that as the value of r increases, the range of stabil- 
ity of the trivial fixed point is the same all throughout, which is — 2 < r < 0. For 
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Time Delay 


Trivial fixed point 


Range of stability 


r = 


Xi = 


-2 < r < 


r = 1 


Xi = (0,0) 


-2 < r < 


r = 2 


Xi = (0,0,0) 


-2 < r < 


T = 3 


Xi = (0,0,0,0) 


-2 < r < 









Table 1 : Time Delay Stability of the trivial fixed point 



the non-trivial fixed point, we can observe that as the r increases, its range of sta- 
bility decreases. From the computations made, it shows numerically the statement 
of Murray that the effect of delay in models is usually to increase the potential for 
instability. Hence, the following proposition is formed from these observations. 



Proposition 3.1. Let 

X„+i = X„ + rX„ (l - ^) (3.1) 

for any r. Then, there are always two fixed points Xi — (0, 0, • ■ ■ , 0), ^2 = 
{K,K,--- , K) and Xi is stable if -2 < r < 0. 



r+1 



r+1 



Proof. Given the dicrete logistic equation with time delay in this form 

^Xn + rX„ (l - ^) (3.2) 
there are always two fixed points, Xi = (0, 0, ■ ■ ■ , 0), X2 = {K, K, - ■ ■ ,K). 

r+1 r+1 

To get the range of stability of the trivial fixed point, we need its Jacobian matrix 
and it is given by 



J{X,) 



/ 1 
10 

■•. 




\ 



\ 









(3.3) 



1 + r / 



(r+1) X (r+1) 
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Since J{X-\) is an upper triangular matrix, then the eigenvalues of it is just its 

diagonal entries, Aj = (2 = 1, 2, ■ ■ ■ , n), At-+i = 1 + r. Applying the stability 
test (cf. [1], etc . . . ), the fixed point Xi = (0, 0, ■ ■ ■ , 0) is stable if -2 < r < 0. 

T+1 

Thus the proof is completed. □ 



4 Stability of the Non-trivial fixed point 

We then use the revised Jury's Conditions for Stability Test in finding the stability 

of the discrete logistic equation with time delay, given by Equation (2.1), at its 
non-trivial fixed point. Having computed for the range of stability of the equation 
when T = 0, 1, 2, 3, we have the following results summarized by Table 2. The 



Time Delay 


Non-trivial fixed point 


Range of stability 


T= 


X2 = K 


< r < 2 


r = 1 


X2 = iK, K) 


< r < 1 


T = 2 


X2 = iK, K, K) 


< r < 0.618034 


r = 3 


X2 = iK,K,K,K) 


< r < 0.445042 









Table 2: Time Delay Non-trivial Fixed Point Range of Stability 

results we have from our computations that is summarized by the table above 
evidently show that as time delay increases, the range of stability of the equation 

determined by the value of r shrinks. According to Murray [2], the effect of 
delay in models is usually to increase the potential for instability. This means 
that as we increase the time delay of the population to respond to changes in the 
environment, the risk of instability increases. With such instability, tremendous 
changes is possible to occur to the population. 

Conjecture 4.1. Let 

= X„ + rX„ (1 - ^) (4. 1) 

for any r. Then there are always two fixed points Xi 



{K, K, - ■ ■ ,K) and X2 is stable if < r < f{T), where 

^ V ' 

T+1 

/(r) : N ^ R is a monotone decreasing function ofr. 
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This conjecture generalizes the behavior of the stability of the discrete logistic 
equation with time delay at its non-trivial fixed point. The proof of the conjecture 
above is given by induction. 

Proof. We begin the proof by assuming that P^^\(A), a polynomial that satis- 
fies /c + 2 Jury Condition. Then afterwards, we will show that the polynomial 
Pn^i^ (A) will be able to satisfy k-\-?> Jury conditions, that is, for any n> k. 
Now, we define Pi%{X) as 



(n+2-fe)-2 



"(71+2-^)- 



(4.2) 



where aS,''^ aj^^^2-fe)-i' «i+2-fe and others are 0. 

^i+i(A) have a value satisfies the /c + 2 Jury Condition given by 



a. 



(k) 

n+2-k 



> % 



(4.3) 



which does not depend on any n> k. Then, Pn+i^\^) is given by 

p(fe+l)/xN Xn+2-(fe+l) , ^(fc+l) X(n+2-(fc+l))-l , _L ^^=+1) X2 
l^J ^ ' H h a(„+2-(fc+l))-2'^ 



"'(n+2-(fc+l))-l^ ^ "n+2-(fc+l) 



, x2 
"(n+l-fe)-2^ 



-I- /)(*^+-^) X -I- /i*-'^"^^'' 



(fe+1) 



(4.4) 



Since the polynomial is the preceding polynomial after k, by the process of 
obtaining the coefficients described in (5.3) or using Table 3 in Appendix, we must 
use the coefficients of the polynomial P^f|.\(A) in order to find the coefficients of 



pi':V)(A). 



n+l 



Then the coefficients of -P^+t^^ W given by as follows: 



(k+i) 



^n+2-k "(n+2-fc)-l 



(fe) 



(k) (k) 

^{n+2-k)-l^0 



'^n+2-k '^(n+2-fc)-2 
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a 



(fe+i) 

(n+l-fe)-2 



(fe) 



a, 







a 



h 

(fc) 

(n+2-fc)-2 



a 



(fc+i) 

(n+l-fc)-l 



(fc) 



a, 







a 



(k) 

(k) 

(n+2-fe)-l 



Ak) 



'^n+2-fc^(n+2-fc)-l 



a, 



(fc+i) 



a. 



(fc) 

ri+2-fc 



a 



(fc) 



(fc) 
(fc) 

''n+2-k 



a, 



(fe) 

n+2- 



(fe)N2 



-kf - 



Due to this, we say that the only coefficients of the polynomial pj:|'_^^^\X) 
which are not equal to zero are Oq'^^^'*, a|^^^^_^^_^, a|j^/2fe- Determining the A; + 3 
Jury Condition, we have 

(4.5) 



ffc+i) 



I'^n+l-fcl > 1^0 



Then, since this condition succeeds k + 2 Jury Condition for any n > k, we 
can say that this condition futher shorten the range of stability of r. And thus, 
proving the conjecture. 

□ 



5 Appendix 

In this appendix, we present Jury Conditions in the following form, which is es- 
sentialy same to Jury Conditons for Stability Test given in [2] (cf . Appendix B in 
[2]). 

Given a Jacobian matrix with dimension (n + 1) x (n + 1), a characteristic 
polynomial of order n + 1 is of the form: 



P(A) = aoA"+^ + aiX" + asA""^ + • • • + a„A + a„+i (5.1) 

where e R, i = 0, 1, • • • ,n + l. 

From this, we have the characteristic polynomial of order n + 1 given by: 



(1) 

n+l 



Pn+iW = 4'^A"+^ + a^h'^ + 4'U"-i + . . . + a^A + a 

where af ^ e R, i = 0, 1, • • • , n + 1, a^^ji 7^ 0. 

The formulas used in defining bk,Ck,--- , ?fe are rewritten as follows: 



(5.2) 
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a 



(2) 



a, 



"it+l 



k = Q,l,2,...,n 



a 



(3) ._ 



(2) (2) 



,(2) 



n— 1 — fc 
(2) 



a 



A; = 0,l,2,...,n-1 



(5.3) 



a 



(n+2) 



a 



a. 



(n+l) ^(n+1) 
(„+l) („+i) 



*0 "fc+l 

Thus, the Jury's stabiUty table is given by: 



A; = 0,l,2 



row 


A° 




A2 


A3 




^„-2 


^n-1 


A" 


^n+1 


1 






"n-l 


"n-2 




4^^ 


4^) 




"o 


2 


"o 




(2) 
^2 








"n-l 




"n+l 


3 


(2) 

an' 


(2) 


(2) 
'^n-2 


(2) 




(2) 
^2 


(2) 

a{ 


(2) 
^0 




4 


(2) 


(2) 


(2) 
^2 


(2) 
«3 




(2) 
^n-2 


(2) 
^n-1 


(2) 
On 
























2n-5 








(n+l) 
^0 












2n-4 


"o 






(n+l) 
^3 












2n-3 


"2 


An+2) 

















Table 3: Revised Jury's Stability Table 
Finally, we rewrite the Jury Conditions for stability given as follows: 

I "n+l I ^ "O 

4. I4?l>"i^^ 

5. |41l>af 

(n + 2). |4"+')|>4"+') 
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Also, we need to define polynomial in which the coefficients are from Table2, 
^i+iO = 2, 3, . . . , n + 2). Using the same notations we had in (2.3), we denote 
the polynomials as follows: 

a?) A" + a?X--' + afr-' + • • • + a^.A + ai^), 
a?) A"- + afh--' + af A-^ + • • • + a£,A + a^l„ 

(5.4) 
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^Si(A) = 



P. 



(n+l) 



n+1 
r,(n+2) 
n+l 



(A) 
(A) 



